An old Method of Jacob! to find Lagrangians 

ob"' M.C. Nucci and P.G.L. Leach* 

o: 

O ' Dipartimento di Matematica e Informatica,Universita di Perugia, 06123 Perugia, Italy, 

e-mail: nucci@unipg.it 



(N 



r- 



00 



c^ 



Abstract 



In a recent paper by Ibragimov [N. H. Ibragimov, Invariant Lagrangians and a new 
method of integration of nonlinear equations, J. Math. Anal. Appl. 304 (2005) 212- 
235] a method was presented in order to find Lagrangians of certain second-order ordinary 

G ' differential equations admitting a two-dimensional Lie symmetry algebra. We present a 

method devised by Jacobi which enables to derive (many) Lagrangians of any second-order 
differential equation. The method is based on the search of the Jacobi Last Multipliers of 
the equations. We exemplify the simplicity and elegance of Jacobi's method by applying it 
to the same two equations as did Ibragimov. We show that the Lagrangians obtained by 

^ , Ibragimov are particular cases of some of the many Lagrangians that can be obtained by 

^T ■ Jacobi's method. 

^' ■ 1 Introduction 

O 

(yQ [ The method of the Jacobi last multiplier [U [H [3l [H [5] provides a means to determine an 

^^ ■ integrating factor of the partial differential equation 

X'- Af = ^a,{xi,...,Xn)-^ = (1) 



or its equivalent associated Lagrange's system 

dxi dx2 dxn 

The multiplier M is given by 

d{f,U;i,L02,---,UJn-l) 
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(9(xi,X2,... ,Xn) 



MAf, (3) 
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and uJi, 



, iOn 



_i are n — 1 solutions of ([T|) or, equivalently, first integrals of ([2]). Jacobi also 



proved that M is a solution of the following linear partial differential equation 
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(5) 



In general a different selection of integrals produces another multiplier, M. An important 
property of the last multiplier is that the ratio, M/M, is a solution of ([1]), equally a first 
integral of ([2]). Indeed, if each component of the vector field of the equation of motion is free of 
the variable associated with that component, ie dui/dxi = 0, the last multiplier is a constant. 
In its original formulation the method of Jacobi last multiplier required almost complete knowl- 
edge of the system, ([U or ([2]), under consideratior(^. Since the existence of a solution/first 
integral is consequent upon the existence of symmetry, an alternative formulation in terms 
of symmetries was provided by Lie [71[8][Kap 15, §5 in the latter]. A clear treatment of the 
formulation in terms of solutions/first integrals and symmetries is given by Bianchi [9j. If we 
know n — 1 symmetries of (ll5p /( fT6]) . say 



i j — y ^ qijyXi, • • • , Xn)Oxj ; 1 — i.,n i, 

Jacobi Last Multiplier is given by M = A^^, provided that A 7^ 0, where 
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There is an obvious corollary to the results of Jacobi mentioned above. In the case that there 
exists a constant multiplier, the determinant A is a first integral. This result is potentially 
very useful in the search for first integrals of systems of ordinary differential equations. In 



^Although we should underline that Jacobi himself found last multiphers for several equations without any 
knowledge of its solutions [H [2 El H [5] . 



particular this feature was put to good use with the Euler-Poinsot system [TOj and the Kepler 

problem |llj . 

The following relationship between the Jacobi Last Multiplier and the Lagrangian [5], [6] 

for a one-degree-of-freedom system 

y" = f{x,y,y'), (9) 

where the prime denotes differentiation with respect to the independent variable x, is perhaps 
not widely known although it is certainly not unknown as can be seen from the bibliography 
in [12]. Given a knowledge of a multiplier, namely a solution of the equation ([5]), i.e. 

A(logM) + |=0. (10) 

then ([8]) gives a simple recipe for the generation of a Lagrangian. The only possible difficulty 
is the performance of the double quadrature. Considering the dual nature of the Jacobi Last 
Multiplier as providing a means to determine both Lagrangians and integrals one is surprised 
that it has not attracted more attention over the more than one and a half centuries since its 
introduction. The bibliography of [12j gives a fair indication of its significant applications in 
the past. In more recent years we have presented the application of Jacobi Last Multiplier to 
many different problems [TOKIIKIlKISlOiaiiailZlIIH]- 

In a recent paper Ibragimov [H] proposed a practical approach to the resolution of the classical 
problem of finding the Lagrangian given a second-order ordinary differential equation. In his 
method Ibragimov introduced the idea of an invariant Lagrangian, and derived Lagrangians 
of two second-order differential equations after lengthy calculations involving integration of 
auxiliary differential equations. For the details of the method the interested reader should 
consult the paper [19] . 

In this paper we exemplify the simplicity and elegance of the forgotten method devised by 
Jacobi for finding Lagrangians by applying it to the same two equations as did IbragimoMj. 
Specifically we obtain Jacobi Last Multipliers, and therefore Lagrangians, of the equations 

(11) 

^Both equations are found in the textbook [2(7. The first is example (f2.27) on page 29f and the second is 
Exercise 12.3 on page 300. 
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which possess the Lie point symmetries 

Ti = 2xdx + ydy, 
and 



r2 = x^dx + xydy 



Si 



vdr — 2d, 



yi 



S2 = xiog(x)a^-2(i + iog(x))a, 



y 



(13) 



(14) 



respectively. We note that both symmetries in (J13p and in (J14p generate a Lie's Type III 
algebra [8j, namely a nonabelian and transitive Lie algebra [9j. 



2 Jacobi Last Multipliers and Lagrangians for (1111) 



The calculation of the Jacobi Last Multiplier requires that the differential equation under 
consideration be written as a system of first-order equations. Thus (fTTIl becomes 
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XUl 



(15) 



with ui = y, and U2 = y' ■ The formula (fTOl) for the last multiplier gives a nonlocal exp[— / u^ dx] 
which is not very useful. However, we do have the route, d?]), through the determinant of the 
vector field and the two symmetries. Thus we have 
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{XUIU2 + X — Ul'^){xU2 — Ul) 



Ul 



(16) 



so that the multiplier is 



M 



Ul 



12 



(17) 



{xUlU2 + X — Ul'^){xU2 — Ul) 

If we integrate M12 twice with respect to U2, then from formula ([8]) we obtain the Lagrangian 



L 



Ul 



12 



-^ {XU2 - Ul) \og{xU2 - Ul) + 



XU1U2 + X — U^ 



X'^ 



log{xUlU2 + X — u1] 



-^ + flix,Ui)u2 + /2(x,Ul), 



(18) 



where /i(x,mi) and f2{x,ui) are arbitrary functions of integration. If we substitute (jlSp into 
the Euler-Lagrangian equation, we obtain the constraint 



dh_dh 
dx dui 



X — u\ 
x:^ui 



(19) 



on the hitherto arbitrary functions /i and /2. This Lagrangian was not found by Ibragimov. 
As it was shown in |16j . |18j . /i,/2 are related to the gauge function g = g{x,ui). In fact, 
we may assume 

f = -^ f = ^ + 2xlog(ni)-^f 

dui ' dx 2x^ 



-ui^ = 


and XU2 - 


- u 


.e.: 








and y' = 


X 



namely the arbitrariness in the Lagrangian (jlSp can be expressed as a total time derivative. 
Such a Lagrangian has been termed 'gauge variant' [21] and is notable in that the presence of 
the arbitrary function g has no effect upon the number of Noether point symmetries [16]. In 
this respect it could be regarded as part of the boundary term in the way Noether put it in her 
formulation of her theorem [23|. The class of Lagrangians described by (jlSp is an equivalence 
class. 

We observe that there are two singularities given by 

XU1U2 + X — ui^ = and xu2 — ui = (21) 

When we solve these two equations, i.e.: 

2/' = -- + -, and y' = ^, (22) 

y X X 

we recover the singular solutions of (jlip associated with the singularities of the Lagrangian 

If we take the Lagrangian ([18]) subject to (fT9]) and calculate its Noether point symmetries, 
we find that there is a single Noether point symmetry which is r2 in (I13p . The corresponding 
integral is 

/= ^ ,. (23) 

XU1U2 + X — uf 

With an integral and a multiplier we can generate a second multiplier by a reversal of the 
property that the quotient of two multipliers is an integral. The multiplier is just 

M, = ^ = '—. (24) 

I XU2 — Ui 

Now we can calculate a second Lagrangian from the multiplier (|24p. and find 

Li = K log(tfi - XU2) H \- fiix,ui)u2 + /2(x,'Ui), (25) 



x^ ■ X 



with the constraint: 

dh_df2^^ ul + x 
dx dui x^uf 

or equally in terms of the gauge function g = g{x,ui) 

dg_ f ^^, x-u l 

OUl OX UlX 



(26) 



/i = ^, /2 = ^ + ^^-?. (27) 
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If we take the Lagrangian Li in (|25p and calculate its Noether point symmetries, we find that 
there is a single Noether point symmetry which is r2 in (J13p . The corresponding integral is / 
in ([23]). 

We can generate many (infinite) different Jacobi Last Multipliers of equation (jlip and 
consequently many (infinite) different Lagrangians. In fact we may take any function of the 
first integral / in (j23p and then its product with either M12 or Mi will generate a new Jacobi 
Last Multiplier. For example, we obtain the following multiplier: 

— 2 XU-iUo + X — u^ , ^ 

M2 = M12-PT = 2^^ ^. 28 

1^ Ui[XU2 — Ui) 

and consequently Lagrangian: 

^Ul-XU2. , . U2{uiU2-2) 

L2 = -2 log(ni -XU2) H h fi[x,ui)u2 + /2(x,ui), (29) 

XUi Ml 

with the constraint: 

^-^ = -A. (30) 

dx dui uf 

or equally in terms of the gauge function g = g{x, ui) 

dg_ f ^ ^ _ i 

OUi OX U 



h = ^. h = -l-^_. (31) 



If we calculate the Noether point symmetries of the Lagrangian L2 in p9|) . we find that 
both Fi and F2 in (jlSp are Noether point symmetries. The corresponding integrals are 

I\ = log ( U\U2 1 9 {u^U2 — XUiUo — 2xuiU2 — x) , (32) 

\ x J uf 

and 

1 fxuiU2 + x-uj\ 

^'=P = ( u, J ' (^^) 



respectively. It is worth noting that both singular solutions obtained in ()22p correspond to these 
integrals taking the particular value of zero, namely, when each integral is a configurational 
invariant [22], [23], we obtain a singular solution. 

One of the two Lagrangians derived by Ibragimov [20j[eq (42), p. 223] for equation (jlip is 
the following 

LNi = —+m- —) log ( ^ - ^1^2) , (34) 

xui Vx^ X / \ X J 

which is a particular case of the Lagrangian Li in ()25p with 

/i = i(-log(ui)+log(x)-l), /2 = ^(log(ni)-log(x)) + — . (35) 

X X^ XUl 



The other Lagrangian of Ibragimov [20] [eq (54), p. 225] is the fohowing 

Ln2 = 2 + ^ ~ 2 h M2 - 2 log -!- - UiU2 , (36) 

MJ X^ X \X Ul J \ X J 

which is a particular case of the Lagrangian L2 in ()29p with 

2 
/i = (nf -rclog(ui) +xlog(x) -x) , 

XUl 

/2 = ^{x^-uj + 2xujlog{ui)-2xujlog{x)). (37) 



3 Jacobi Last Multipliers and Lagrangians for (1121) 

The system of first-order differential equations corresponding to (I12p is 

u'l = U2 

4 = -— + e"\ (38) 

X 

In this case the application of formula ([5|) or equivalently ([TO]) does produce a multiplier. It is 

Mo = x (39) 

from which we obtain the Lagrangian 

2 
Lo = ^- + fiix,ui)u2 + f2{x,Ui) (40) 

with the constraint on the two functions of integration being 

^-^ = -xe'^S (41) 

dx dui 

or equally in terms of the gauge function g = g{x, ui) 

h = ^, /. = |^ + xe-. (42) 

GUI OX 

This Lagrangian admits one Noether's symmetry namely Si in (|14p and yields the following 
first integral: 

/o = 4xu2 + U2^x'^ - 2e"ix2. (43) 



We use the two symmetries Si, S2 in (HH) and the vector field of the system ([38]) to obtain 
a second multipher. The matrix is 



Mat 



12 
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U2 


-^+e- 


xlog(2;) 


-2(l + log(x)) 


2 
n2(l + log(x)) 

Jb 


X 


-2 


-U2 



and the corresponding multiplier is 

Mi2 = 



Thus formula ( 
L12 



4 + 4:XU2 + it2^a^^ — 2e"ix-^ 
yields the following Lagrangian 

- log {-XU2 - 2 - \/2xe"i/2') - — log ' ^^^ 



+ 2 + ^/2xe"i/2~ 



2x 



+ 2 - ^/2xe"i/2 



XM2 



^/2(xn2 + 2)^ /xn2 + 2 + \/2xe"i/M 
+ — . n .. /n log r- „ ,^ + hu2 + /2 



with the constraint 



4a,2g«i/2 -^ I a;^2 + 2 - \/2xe"i/2 
dfi 9/2 



Ox 9n 
or equally in terms of the gauge function g = g{x, ui) 



0. 



r dg 
h — 7^ — , 

OUl 



■/'2 — ■^■ 
ox 



(44) 



(45) 



(46) 



(47) 



(48) 



In a curious repetition of the situation with (jlip we find that the two Lagrangians Lq in 
(j40p and L12 in (j46p have the same Noether point symmetry Ei in (J14p and lead to what is 
functionally the same integral Iq in (j43p . Both Lagrangians were not found by Ibragimov. 
Indeed Ibragimov did not look for Lagrangians of equation ()12p admitting one Noether point 
symmetry. 

Since we have two multipliers, we can obtained a first integral given by their ratio, namely 



M12 



Ml ^2 



-(xM2 + 2)^ + 2e"ix 



-In - 4. 



(49) 



We can generate many (infinite) different Jacobi Last Multipliers of equation (J12p and 
consequently many (infinite) different Lagrangians. In fact we may take any function of the 



first integral /q in (j43p and then its product with either M12 or Mq will generate a new Jacobi 
Last Multiplier. For example, we obtain the following multiplier: 

M, = - ^^° =- . " (50) 

and consequently the following Lagrangian: 



, 2\ / J{xu2 + 2)2 - 2e"ia;2 + a;«2 + 2 

-^2 = - n2 + - log ' '' 



+ v/(a;n2 + 2)2-2e«ix2 + /i(x,ui)u2 + /2(x,ni), (51) 

with either the constraint (I47p or (1481) . This Lagrangian admits both Si and S2 in (J14p as 
Noether point symmetries, and the corresponding integrals are 

h = V/o + 4 = V(a;w2 + 2)2-2e"ix2, (52) 

and 



/7 T^5 7 . N / \/(xU2 + 2 2 - 2e"ix2 + OT2 + 2 \ , , 

h = x/(a;w2 + 2)2-2e"ix2log(x) + 21og ^^^^ '-^ , (53) 

\ V2xe"i/"' / 

respectively. 

The Lagrangian of Ibragimov [20] [eq (90)j_p. 234] is a particular case of the Lagrangian L2 
in (j5ip with the gauge function equal to zerco. 



The last multiplier M12 in (j45p becomes singular if 

2 



y' = --±V2ey/'^. (54) 



Equation (fT2]) is satisfied by each of the first-order equations in (fMl) and so we obtain the two 
singular solutions 

y = X fc =F V2x) (55) 



thereby supplementing the results given in [19j. 
^In [20] there are some missprints. 



4 Final Remarks 

When one seeks a Lagrangian of an elementary equation, it is usually possible to guess the 
form of at least one Lagrangian. In the case of not so elementary equations an approach 
using guesswork is likely to lead to frustration. Consequently any development which can 
replace guesswork or intuition by a well-defined procedure is to be welcomed. Usually there 
is a price to pay for the elimination of guessworlQ. In this paper we have considered two test 
equations proposed by Ibragimov to illustrate his concept of the use of invariant Lagrangians to 
provide a new method for the integration of nonlinear equations. We have demonstrated that 
some quite old knowledge is available for a successful resolution of the same problems. The 
combination of the concept introduced by Jacobi in his last multiplier, the application by Lie 
of his ideas of invariance under the transformations generated by continuous groups and the 
specialisation to the Action Integral by Noether provides us with a very powerful and simple 
tool for the resolution of ordinary differential equations which possess a reasonable amount 
of symmetry. We have seen in the two examples considered here that they provide richer 
results when considered from a more classical viewpoint. The Jacobi last multiplier yields 
more general Lagrangians than those find by Ibragimov, and many more can be generated. 
One could consider that the combination of Jacobi and Lie gives sufficient material to deal 
with these equations. In that sense it could be argued that the theorem of Noether is already 
implicit in the work of Jacobi and Lie. However, we did see that further results were available 
to us by an application of Noether's Theorem to the information already obtained. 
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